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SUMMARY 


This  report  describes  the  continuation  of  a  work  programme  [1]  for  the  design,  manufacture  and  testing  of 
a  passive  load  gust  alleviation  device  for  use  on  a  sensorcraft  structure.  The  testing  of  a  prototype  wind 
tunnel  model  is  described  along  with  the  development  of  an  optimised  aeroelastic  scaling  approach  which 
will  be  used  for  the  design  of  a  future  wind  tunnel  model. 


1.  INTRODUCTION 

Increasing  the  aspect  ratio  of  an  airplane  reduces  the  induced  drag  and  therefore  also  fuel  consumption, 
which  is  a  great  importance  for  aerospace  vehicles  such  as  Sensorcraft  that  are  designed  for  long 
endurance  missions.  Unfortunately  a  high  aspect  ratio  also  increases  the  wing  bending  moment  and  in 
case  of  the  Sensorcraft  reduces  the  critical  buckling  speed.  The  design  root  bending  moment  stems  either 
from  manoeuvres  or  gusts  and  this  gives  the  critical  loading  case  and  consequently  leads  to  higher 
structural  mass  to  withstand  the  loads. 

One  way  to  reduce  these  loads  is  by  means  of  active  manoeuvre  load  control  or  gust  alleviation.  This  is 
performed  by  the  movement  of  the  control  surfaces  through  the  flight  control  system  but  requires  a  highly 
redundant  system  and  associated  controls  technology  and  power  requirements. 

The  alternative  approach  is  to  implement  a  passive  load  alleviation  system  which  does  not  require  a 
control  system  and  can  be  made  failsafe  without  redundancy  requirements.  The  earlier  part  of  this  project 
investigated  the  use  of  a  passive  device  and  demonstrated  its  feasibility  analytically  on  a  baseline 
Sensorcraft  model.  It  was  found  that  it  would  be  possible  to  reduce  the  design  loads  and  overall  structural 
weight  if  such  a  device  was  used. 

This  report  describes  the  work  continuing  the  development  of  a  passive  load  alleviation  system  from  July 
2006  to  June  2006  under  EOARD  contract  FA8655-05- 1-3006. 

A  simplified  1:40  elastic  Sensorcraft  model  wing  with  passive  gust  alleviation  device  was  designed, 
manufactured  and  tested  in  the  low  subsonic  wind  tunnel  of  MACE  (Manchester  University).  The  wind 
tunnel  tests  show  a  reduction  of  wing  deflection  up  to  about  50%  by  using  the  gust  alleviation  device. 
Further  effort  was  devoted  towards  developing  an  aeroelastic  scaling  optimisation  process  for  the  design 
of  the  wind  tunnel  model  that  is  to  be  tested  in  the  next  part  of  the  project. 
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Figure  1 .  Planform  Arrangement  of  Sensorcraft  and  Baseline  Configuration  Parameters[2] 


2.  BASELINE  SENSORCRAFT  STRUCTURE 

The  Sensorcraft  planform  and  dimensions  used  to  provide  the  baseline  structure  for  the  scaled  prototype 
wind  tunnel  model  are  presented  in  Figure  1  [2],  It  was  decided  to  apply  a  geometric  scaling  factor  of  40 
in  order  to  fit  into  the  wind  tunnel  that  was  to  be  used  for  testing  of  the  prototype. 


2.1.  MODEL  DIMENSIONS  FOR  1:40  MODEL 


As  can  be  seen  in  Figure  1,  the  total  40th  scale  model  has  a  span  of  —  =  0.85  m  and  the  outer  wing  gust 

g 

device  —  =  0.2  m. 

40 

2.2.  ATTACHMENT  SPRING  FOR  PASSIVE  GUST  DEVICE 

A  schematic  representation  of  the  gust  device  can  be  found  is  shown  in  Figure  2.  Essentially  the  wing  is 
‘cut’  towards  the  outer  tip  and  both  parts  are  joined  via  a  torque  tube  along  the  leading  edge.  As  the  outer 
part  is  attached  well  forward  of  the  aerodynamic  centre  on  the  quarter  chord,  the  effect  of  any  vertical  gust 
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is  to  produce  a  nose  down  motion  in  the  wing  tip,  thereby  counteracting  the  effect  of  gust.  As  the  gust 
device  is  at  the  wing  tip,  the  load  alleviation  will  have  a  significant  effect  on  the  resulting  loads. 


Out  board  section 


Figure  2.  Schematic  of  Gust  Load  Alleviation  Device 


2.2.1 .  Initial  Modeling  Approach 

The  initial  approach  to  confirm  the  underlying  idea  behind  the  passive  device  was  the  application  of  a 
discrete  point  loading  on  the  baseline  FE  wing  model  which  was  modified  by  the  introduction  of  a  ‘cut’  to 
include  the  passive  device.  This  single  point  of  action  on  the  outer  wing  was  defined  to  represent  the 
application  of  a  gust.  The  resulting  washout  angle  was  determined  and  the  consequent  discrete  gust  load 
after  washout  could  then  be  defined. 


2.2.2.  Definition  of  a  Discrete  Gust  and  Point  of  Action  on  the  Outer  Wing 

With  a  gust  of  60  ft/sec  (=18.15  m/s)  at  65  m/s  S.L.  we  get  a  gust  angle  of  0.279  rad  or  16  degrees. 
Taking  the  gust  force  as: 

F7  =  Caa  a  / 57.3  q  F 


leads  to  a  pressure  of 


6521.225 

2 


=  2.6  KPa  and  hence  the  force  becomes 


4.93  x  0.279  x  2.6  x  20  =  71 .5  KN 
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2.2.3.  Modification  of  Finite  Element  Model 


The  outer  wing  was  only  attached  at  the  leading  edge  as  shown  in  Figure  3  along  with  the  position  of  the 
discrete  gust  force.  The  modification  was  only  applied  on  the  two  ribs  and  the  ‘nose  spar’.  The  resulting 
washout  angle  was  calculated  from  the  z  displacements  of  nodes  610  and  591  and  648  and  629  as  shown 
in  Figure  4.  The  resulting  washout  angle  was  found  to  be  -6.25  degrees  or  0.109  rad.  This  result  can  be 
altered  by  changing  the  torsional  stiffness  of  the  front  spar.  A  more  flexible  spar  leads  to  a  greater 
washout  and  vice  versa. 


Figure  3.  FE  Model  Showing  “Cut”  and  Point  of  Discrete  Loading 


Figure  4.  Stiffened  structural  elements,  “leading  edge”  elements  t=40mm,  ribs  elements  t=5mm 


From  this  result  the  effective  torsional  stiffness  can  be  deduced  using  an  arm  of  0.7  m  (distance  of  outer 
wing  pivot  to  quarter  chord  axis)  such  that 


K  =  ■ 


M 


71.5x0.7 

0.109 


=  459.2KNm/rad 
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The  equivalent  gust  angle  for  a  60  ft/sec  gust  =18.15  m/sec  is  derived  as 


18.15 

65 


=  0.279  rad  =  16  degrees  at  65  m/sec. 


Therefore  an  alleviation  factor  of 


0. 108eiastjc  — 

0-279rigid 


can  be  calculated. 

For  the  design  of  the  wing  tunnel  model,  we  want  to  achieve  a  0.279  rads  (16  degrees)  rigid  gust  angle  at 
25  m/sec,  which  gives 


252x  1.225 
2 


0.38KPa 


and  a  vertical  force  of 


4.93  x  0.279  x  0.38  x  20  =  1 0.45  KN  (same  wing  @  lower  speed) 


To  get  0.109  rad  for  the  elastic  sensorcraft  at  25  m/sec,  the  torsional  stiffness  must  be 


459.2x10.45 

TL5 


=  67.1  KNm/rad 


2.2.4.  Wind  Tunnel  Model  1 :40  Scale  Spring 


Taking  67.1  KNm/rad  and  multiplying  it  with 


7  i  V 


v40y 


of  Table  1  a  stiffness  of 


67.1  x  0.025  =  1 .049  Nm/rad  results. 
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Dimensions 

Scaling  Factor 

Units 

Length 

0.025 

Mm 

Stiffness 

0.025 

N/m 

Torsional  Stiffness 

0.0253 

Nm/rad 

Frequency 

40 

Hz 

Velocity 

1.0 

m/s 

TABLE  1  SCALE  FACTORS 


3.  MODEL  DESIGN 

In  Figure  5,  a  picture  of  wing  device  /  full  wing  system  is  shown  including  rib  labeling. 

Connector 


I  t 


t 


1 1  i 


Soot  Rib  1  Rib  2  Rib  3  Rib4  Bib  5 

Figure  5  Illustration  of  Wing-Device/Full  Wing  system  with  rib  labelling 
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The  dimensions  of  the  wing  are  presented  in  Figure  6. 


3.1.  MANUFACTURE  OF  WIND  TUNNEL  WING 

3.1 .1 .  Manufacture  of  the  Main  Components 

It  was  decided  that  for  the  length  of  the  %  wing  and  device,  the  spars  would  run  straight  through  each  rib 
from  one  end  to  another.  This  would  help  the  model  to  maintain  its  structure  when  it  was  not  subjected  to 
any  loads.  In  order  generate  the  slots  through  the  ribs,  a  hole  was  initially  drilled  at  either  end  of  each  slot 
and  was  followed  by  the  removal  of  the  remaining  material  in  between  using  a  CNC  miller,  as  shown  in 
Figures  7  and  8. 
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The  ribs  were  machined  from  a  rectangular  aluminum  bar  using  the  same  CNC  miller  as  above  so  that 
they  had  an  aerofoil  profile.  The  profile  chosen  was  the  NACA  0012  shown  in  figure  9.  The  code  ‘0012’ 
indicates  that  aerofoil  has  no  camber  and  a  maximum  thickness  equal  to  12%  of  its  chord. 
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In  milling  the  profile  the  cutter  traced  the  profile  ten  times  so  as  to  cut  tirough  the  15  mm  thick  block 
gradually,  i.e.  with  each  cycle  the  cutter  increased  its  depth  by  1 .5  mm  as  shown  in  figures  10  and  1 1 . 


Material  for  the  aluminum  rectangular  spars  were  ordered  in  a  long  bar  so  that  it  could  be  cut  down  to  the 
correct  length. 
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Holes  for  the  connector  were  then  drilled  in  the  LE  and  TE  spars  of  the  wing  and  device  on  the  adjacent 
end  faces  and  the  connectors  were  cut  to  length  from  a  silver  steel  rod.  Following  analytical  work,  the 
desired  diameter  of  the  connector  was  known.  However,  due  to  material  limitations,  it  was  not  possible  to 
obtain  rods  of  the  required  size.  The  closest  two  rods  were  found,  one  slightly  larger  and  the  other 
smaller,  from  which  the  smaller  diameter  rod  was  chosen  for  the  current  model.  By  choosing  the  smaller 
diameter  rod  the  stiffness  of  the  connector  is  less  than  that  of  the  desired  diameter. 

Prior  to  assembly  of  the  model,  holes  were  drilled  and  tapped  in  four  positions  along  the  top  surface  of  the 
ribs,  corresponding  to  where  the  spars  passed  through.  Grub  screws  would  then  be  used  as  an  extra 
precaution  to  avoid  unwanted  movement  of  the  spars  and  ribs  relative  to  each  other.  In  addition,  a  further 
two  holes  were  drilled  on  the  bottom  surface  of  the  third  and  fourth  ribs,  as  well  as  the  LE  and  TE  spars, 
corresponding  to  where  the  connector  passes  through.  Similarly,  grub  screws  were  used,  but  with  greater 
significance,  to  fix  the  connector  and  hence  prevent  movement  at  those  points. 


3.1 .2.  Manufacture  of  the  Balsa  Wood  Leading  Edge 

As  the  spars  are  rectangular,  an  aerofoil  shaped  leading  edge  was  made  from  balsa  wood.  This  is  to 
provide  support  for  the  solarfilm  to  be  later  wrapped  around  the  wing  and  device,  which  then  constitutes 
the  skin  of  the  model,  and  to  ensure  that  the  air  passes  over  the  model  correctly.  The  balsa  wood  leading 
edge  was  made  by  gluing  rectangular  strips  of  balsa  wood  to  the  spars  and  filing  down  to  shape,  using  the 
ribs  as  a  guide,  as  shown  in  Figure  12. 
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3.1 .3.  Enabling  Fixture  of  Model  to  Wind  T unnel 

The  model  is  attached  to  the  wind  tunnel  via  a  plate  interface  shown  in  Figure  13.  The  plate  is  fixed  to  the 
wind  tunnel  by  two  screws  but  can  be  adjusted  so  that  it  is  angled  at  certain  degrees.  Three  blind  holes 
were  then  drilled  in  the  ‘root’  rib,  one  for  a  threaded  rod  and  the  remaining  two  were  for  locating  pins. 
These  pins  were  positioned  forward  and  aft  of  the  threaded  rod,  which  was  positioned  in  between  the  two 
inner  spars.  The  pins  located  into  the  blind  holes  in  the  plate  to  ensure  that  the  wing  is  always  correctly 
fitted  in  the  same  position  on  the  plate.  By  adjusting  the  angle  of  the  plate,  the  wing’s  root  angle  of  attack 
is  also  adjusted.  Finally,  the  threaded  rod  fitted  through  the  plate  and  wind  tunnel  and  could  then  be  bolted 
on  to  the  wind  tunnel  wall. 


FIGURE  13.  PLATE  INTERFACE  INCLINES  AT  2s 


3.1 .4.  Skinning  of  the  Model 

Once  the  model  had  been  modified  so  that  it  was  attachable  to  the  wind  tunnel,  the  %  wing  was  re¬ 
assembled  and  the  connectors  removed  from  the  %  wing  and  device.  Thin  sheets  of  balsa  wood  were  cut 
to  tightly  fit  in  between  each  rib  at  the  trailing  edge  and  were  also  held  in  place  by  glue.  The  purpose  of 
these  sheets  was  to  prevent  sticking  of  the  solarfilm  to  itself  from  the  top  and  bottom  surfaces  of  the 
model,  as  this  would  not  yield  an  acceptable  aerofoil  wing. 

The  device  was  then  wrapped  in  solarfilm  and  an  iron,  with  gradually  increasing  heat,  was  passed  over 
the  surfaces.  As  the  solarfilm  heated,  it  contracted  and  shrank  to  fit  the  shape  of  the  device.  The  same 
procedure  was  followed  for  the  %  wing,  however,  wrapping  of  this  wing  was  much  more  difficult  because 
the  rod  and  pins  were  not  removable.  It  was  important  to  gradually  increase  the  heat  and  ensure  that  the 
iron  was  not  on  the  film  for  any  length  of  time,  or  that  the  film  would  not  overly  stretch,  otherwise  damage 
to  the  film  was  likely  to  occur.  This  could  result  in  ripping  of  the  film  and  the  process  would  have  to  be 
repeated.  The  completed  wing  with  gust  alleviation  device  is  shown  in  Figure  14. 
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Following  the  skinning  of  the  model,  it  was  then  necessary  to  make  small  holes  in  the  solarfilm  at  locations 
where  the  grub  screw  required  accessing  in  order  to  fix  or  release  the  connector.  Prior  to  making  these 
holes  tape  was  attached  to  the  surface  of  the  model  at  the  relevant  locations  to  prevent  the  holes  made  in 
the  film  from  spreading. 


4.  WIND  TUNNEL  TESTS 
4.1.  Connector  Stiffness 


The  wing  and  the  gust  device  are  connected  by  a  cylinder  which  is  fixed  to  both  parts  by  grub  screws.  In 
order  to  get  the  desired  stiffness  the  cylinder  diameter  is  only  3mm  and  its  length  was  19mm. The  material 
chosen  was  aluminium. 


The  stiffness  can  be  calculated  for  a  cylinder  as: 


G.Jt  _GA^_ 

1  1x32 


G  Aluminium  =  80000 


N 


mm 


Therefore  for  a  cylinder  of  0.019  m  and  0.003m  diameter  we  get 


K§  = 


28  x  109  x  0.0034  X7t 
0.019x32 


=  1 1 .7  Nm/rad 


The  problem  with  such  a  design  is  that  the  diameter  enters  the  equation  with  the  fourth  power.  The  grub 
screws  which  fix  the  connector  intrude  in  the  cylinder  and  therefore  reduce  the  diameter. 
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Assuming  a  diameter  of  0.0015m  to  allow  for  the  attachment  we  get  0.73  Nm/rad  which  is  very  close  to 
the  1 .05  Nm/rad  which  would  give  approximately  60%  load  alleviation  due  to  the  elastic  twist. 

4.2.  Static  Wind  Tunnel  Tests 

The  tests  were  performed  with  a  root  incidence  of  3  degrees.  In  hindsight  it  might  have  been  better  to 
choose  a  larger  value  due  to  the  presence  of  structural  non-linearities  that  occur  do  to  the  small 
deflections  that  occurred.  The  wing  was  tested  both  with  device  operational  and  also  a  baseline  case 
where  the  trailing  edge  was  fixed  to  the  wing,  eliminating  any  effects  of  the  gust  alleviation  device.  Initially 
the  wing  was  tested  for  varying  tunnel  speeds  and  static  deflections  measured. 

In  figure  15,  deflections  of  rib  3  are  plotted  for  the  whole  wing. 


Figure  16  shows  deflections  at  rib  3  with  the  gust  alleviation  device  working 
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In  Table  2,  the  rib  3  trailing  edge  deflections  are  summarised.  There  are  non-linearities  present  which 
mainly  stem  from  fixing  the  connector  to  the  wing  (grub  screws).  Overall  a  reduction  of  deflection  is 
shown  when  the  gust  device  is  used.  For  the  highest  speed  the  reduction  in  wing  deflection  is  about  53% 
which  correlates  well  with  the  theoretical  calculations  of  the  scaled  1 :40  wind  tunnel  model. 


Full  Wing 

Wing  Device 

%  Difference 

0 

0.1 

5 

0.22 

10 

0.46 

0.27 

42 

15 

0.88 

20 

1.28 

0.92 

27 

25 

3.15 

1.47 

53 

TABLE  2  RIB  3  TE  DEFLECTIONS  (mm)  for  Various  Speeds 


4.3.  Dynamic  Wind  Tunnel  Tests 

A  set  of  gust  vanes  were  manufactured  which  are  able  to  produce  computer  controlled  vertical  gusts  into 
the  tunnel  -  sinusoidal,  “1-cosine”  and  random  gusts  can  be  produced.  Figure  17  shows  the  gust  vanes 
in  the  tunnel  and  the  model. 
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Initial  investigations  concentrated  upon  the  use  of  sinusoidal  gust  excitation,  it  was  possible  to  tune  the 
gust  frequency  in  order  to  excite  the  first  bending  mode  and  to  achieve  significant  deflections.  Figures  18 
and  19  shows  the  wing  tip  deflections  with  the  gust  device  off  and  on.  It  can  be  seen  that  when  the  device 
is  off,  the  leading  and  trailing  edge  deflections  are  very  similar,  indicating  very  little  twist.  With  the  on,  the 
trailing  edge  deflects  more  than  the  leading  edge,  which  shows  that  it  works,  and  the  deflections  of  the 
leading  edge  are  less  indicating  that  the  bending  motion  is  reduced. 


Figure  18.  Leading  and  Trailing  Edge  Deflections  Device  Off 
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m  19mb  Doric*  ON 


Tim*  (a) 


Figure  19.  Leading  and  Trailing  Edge  Deflections  Device  On 


5.  Design  of  Wind  Tunnel  Sensorcraft  Model 

The  next  part  of  this  work  involves  the  design,  manufacture  and  testing  of  a  larger  half  Sensorcraft  model 
to  be  tested  in  the  9’x7’  tunnel  at  the  University  of  Manchester.  This  was  initially  going  to  take  the  form  of 
a  1 :20  model  however  following  discussions  with  Dr  Max  Blair  and  co-workers  the  baseline  structure  has 
changed,  and  will  now  follow  the  Boeing  design  shown  in  Figure  20.  The  gust  alleviation  device  will  be 
incorporated  into  the  structure,  and  in  the  same  way  as  the  prototype  wing  discussed  above,  tests  will  be 
performed  with  and  without  the  device  being  activated. 

It  is  planned  that  a  PhD  student  from  AFIT  will  spend  some  time  in  Manchester  assisting  with  the  tests, 
also  investigating  the  non-linear  geometric  behaviour  at  high  static  loading  cases. 

There  were  a  number  of  difficulties  in  actually  defining  the  parameters  of  the  model  (we  were  not  allowed 
to  have  access  to  the  FE  model)  and  consequently  only  some  initial  analysis  has  been  performed. 
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Figure  20.  Baseline  Structure  for  Future  Work 


6.  Aeroelastic  Scaling 

6.1.  Objectives 

•  To  define  a  dynamically  scaled  aeroelastic  model  that  responds  identically  to  a  full-scale  design 
with  respect  to  chosen  scale  factors  (e.g.  geometric  scaling,  velocity  scaling  etc).  This  goal  can  be 
achieved  on  the  basis  of  a  set  of  non-dimensional  aeroelastic  equations. 

•  To  eliminate  the  reliance  of  the  BAFI  [3]  scaling  approach  on  equivalent  beams 

6.2.  Approach 

•  Based  on  wind-tunnel  parameters,  identify  maximum  geometry  scale  and  maximum  velocity  scale. 
For  a  dynamically  scaled  test,  will  also  need  to  identify  a  third  scale  (e.g.  mass); 

•  Identify  the  governing  non-dimensional  equations  of  motion  (EOM)  that  will  be  numerically 
identical  for  both  the  scaled  and  full-scale  designs.  The  choice  of  degrees  of  freedom  must  serve 
both  models.  The  common  degrees  of  freedom  may  be  reflected  on  a  global  scale  with  modal 
coordinates  or  a  local  scale  with  influence  coefficients  interpolated  at  common  points. 

•  Start  the  scaled  model  design  process  by  defining  the  flexible  structure.  This  layout  could  be  a 
lattice  of  beams  (but  not  necessarily  -  could  be  ribs  and  spars,  or  perhaps  an  equivalent  plate). 

•  Size  the  scaled  model  structure  according  to  the  scaling  parameters.  This  step  depends  on  the 
form  of  the  governing  aeroelastic  equations  of  motion. 

6.3.  Non-Dimensionalised  Aeroelastic  Equations  of  Motion 

In  order  to  achieve  a  non-dimensional  set  of  equations,  we  need  to  select  3  base  units  of  measurement 
for  dynamic  analysis,  typically  mass,  length  and  velocity  (length  /  time)  although  force,  length  and  velocity 
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could  also  be  chosen.  For  static  analysis  we  only  need  2  base  units  of  measurement  -  force  and  length 
(the  time  dependent  term  is  removed). 

The  following  equations  assume  we  have  reduced  the  equations  of  motion  to  a  minimum  of  two  degrees 
of  freedom  that  represent  displacement  (length)  and  rotation  (non-dimensional).  These  two  degrees  d 
freedom  are  consistent  with  most  commercial  finite  element  codes. 

Taking  a  coupled  FE  /  Aerodynamic  equations  of  motion  (ignoring  structural  damping)  such  that 
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Vector  of  translational  degrees  of  freedom 
Vector  of  rotational  degrees  of  freedom 
Block  matrix  terms  in  inertia  /  mass  matrix 

Block  matrix  terms  in  stiffness  matrix 

Reference  length  (semi-chord) 

Block  matrix  Aerodynamic  terms  (complex  for  unsteady  motions) 


In  terms  of  dimensions  of  the  fundamental  quantities  (M,  L,  T)  equation  (1)  has  the  form 
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Thus,  according  to  Newtonian  mechanics,  we  have  three  base  scaling  parameters  of  mass,  length  and 
time. 


6.4.  Base  Scaling  Parameters 


Geometric  Scaling:  Defining  the  non-dimensional  coordinate  system  (£,<|))  such  that 

x  =  be  and  0  =  4>  (3) 
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Time  Scaling:  non-dimensional  time  x  =  — .  Thus,  the  base  time  unit  is  indirectly  scaled  in  terms  of  a 

b 

(surrogate)  velocity  and  geometric  scale  -  which  reflects  the  requirement  that  the  scaled  model  operate 
within  the  wind  tunnel  operational  limits.  Use  the  chain  rule  to  non-dimensionalize  time  derivatives  that 
are  present  in  the  governing  equations  of  motion. 


V2 

x  =  x  =  —  ? 

b 


(4) 


where  (.)  denotes  differentiation  w.r.t  time  and  (‘)  denotes  differentiation  w.r.t  non-dimensional  time. 


Mass  Scaling:  There  are  a  number  of  options  for  scaling  the  mass.  These  are  discussed  below.  Mass 
can  be  scaled  directly  (e.g.  total  mass  of  the  airplane).  Mass  can  also  be  scaled  in  terms  of  a  surrogate 
density  scale.  Also,  since  Force  is  mass  times  acceleration,  units  of  force  can  be  used  as  a  derivative  (or 
surrogate)  scaling  parameter  for  mass. 


Equation  (2)  will  be  rendered  non-dimensional  according  to  the  three  scaling  parameters. 


6.5.  Mass  Scaling 

The  first  approach  will  consider  using  some  reference  mass  mr  as  part  of  the  scaling  process. 

Non-Dimensionalised  Inertial  Force 

Defining  non-dimensional  inertia  terms  (denoted  by  the  overbar) 

Mu  =mrMu  M12=mrbM12  M21=mrbM21  M22=mrb2M22  (5) 

where  mr  is  some  reference  mass  (e.g.  total  mass  of  the  full  structure)  then  the  force  and  moment 
resulting  from  the  mass  and  inertial  terms  can  be  written  as 
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where  F  and  M  are  non-dimensional  forces  and  moments  respectively. 

Non-Dimensionalised  Stiffness  Force 

The  process  for  the  stiffness  terms  follows  the  same  procedure  as  that  shown  above  for  the  inertia  terms. 
Defining 


mi  \  -  _  mi  V  - mi  v7  '  _  _ 

K„=-^KU  K12=^-K12  K21=^-K21  K22=mrV2K22(8) 

then  the  force  and  moment  resulting  from  the  stiffness  terms  can  be  written  as 
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Non-Dimensionalised  Aeroelastic  Equations 

Introducing  transformations  (3)  and  (4)  into  equation  (1)  and  pre-multiplying  by 
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gives  the  redefined  system  equations  as 


(11) 
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Note  that  the  Q  matrix  now  becomes  fully  non-dimensional. 


6.6.  Density  Scaling 

A  second  approach  will  consider  using  the  air  density  as  part  of  the  scaling  process 

Non-Dimensionalised  Inertial  Force 

Defining  non-dimensional  inertia  terms  (denoted  by  the  overbar) 

Mn=pb3Mn  M12  =  pb4M12  M21=pb4M21  M22  =  pb5M22  (13) 

where  p  is  the  air  density  and  b  is  some  defined  length,  then  the  force  and  moment  resulting  from  the 
mass  and  inertial  terms  can  be  written  as 
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Non-Dimensionalised  Aeroelastic  Equations 

Introducing  transformations  (13)  and  (16)  and  pre-multiplying  equation  (1)  by 
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gives  the  redefined  system  equations  as 
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Again,  the  Q  matrix  now  becomes  fully  non-dimensional.  It  is  arguable  that  this  is  a  better  form  to  use  as 
there  is  only  a  length  scaling  that  needs  to  be  changed,  however,  there  are  now  no  longer  any  density  or 
mass  terms  directly  visible  in  equation  (20). 


6.7.  Scaling  of  the  Aeroelastic  Equations 

We  have  arrived  at  two  different  non-dimensional  sets  of  aeroelastic  equations  (12)  and  (20).  These 
equations  can  be  used  to  determine  the  system  behaviour  for  any  combination  of  fundamental  quantities 

p,b,mr  or  p,b,V  respectively  using  the  scaling  relationships  (5)  and  (8)  or  (13)  and  (16)  respectively. 
If  the  full  size  and  model  structure  (physical  construction  and  FE  /  aerodynamic  models)  are  exactly  the 
same,  then  these  scaling  relationships  determine  the  scaled  model  characteristics. 

However,  in  practice  the  model  structure  is  not  constructed  in  the  same  way  as  the  full  size  aircraft,  and 
the  finite  element  and  aerodynamic  models  will  also  not  be  the  same  therefore  the  above  scaling 
relationships  can  be  used  to  provide  “target”  characteristics  for  the  scaled  model.  For  the  rest  of  this 
section  the  mass  scaling  approach  will  be  used  but  the  density  scaling  approach,  or  any  other,  could  be 
used  as  well. 


Wind-Off  Natural  Frequencies  and  Mode  Shapes 

The  natural  frequencies  and  corresponding  mode  shapes  of  the  full  scale  system  are  obtained  from  the 
eigensolution  of  the  matrix 


D<I>  =  <I>n  (21) 

where 
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and  the  natural  frequencies  have  units  of  rad  s'1  whereas  the  mode  shapes  are  non-dimensional. 

In  terms  of  the  non-dimensionalised  equations  of  motion  (taking  the  mass  scaling  route),  the  natural 
frequencies  and  corresponding  mode  shapes  are  obtained  from  the  eigensolution  of  the  matrix 
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and  the  natural  frequencies  now  have  units  of  rad  x_1 .  Unfortunately  at  zero  airspeed  the  scaling  has  no 
meaning  and  therefore  the  dynamic  analysis  must  be  performed  with  the  inclusion  of  the  velocity 
dependent  aerodynamic  terms. 


Wind-on  Natural  Frequencies  and  Mode  Shapes 

Considering  the  system  when  there  is  some  airspeed,  the  non-dimensionalised  frequencies,  dampings 
and  mode  shapes  are  found  from  the  solution  of  equation  (25). 
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Although  the  structural  terms  remain  the  same  whatever  the  scaling,  the  aerodynamic  terms  do  change 
depending  upon  the  triple  (rr>,  ?,  b).  Remember  that  the  frequencies  that  are  obtained  are  in  non- 

dimensional  time  and  thus  are  scaled  by  b/y  and  that  the  mode  shapes  are  complex. 

Static  Aeroelastic  Deflections 

The  static  aeroelastic  deflections  depend  upon  there  being  some  initial  incidence  on  some  parts  of  the 
lifting  surfaces.  In  terms  of  the  non-dimensionalised  equations  (25)  we  get 
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where  <|)0  are  a  set  of  initial  angles  of  incidence.  The  deflections  (;,<()  that  are  found  are  non- 
dimensionalised. 


6.8.  Practical  Application  of  Aeroelastic  Scaling 

In  practice,  the  problem  is  defined  as  to  find  an  equivalent  aeroelastically  scaled  model  structure  given  the 
FE  and  aerodynamic  models  of  a  full  scale  structure  at  specified  flight  conditions.  Although  the  geometric 
scaling  holds  for  the  planform  of  the  model,  the  internal  structure  is  likely  to  be  different,  possibly  not  even 
made  of  the  same  material  or  construction  technique;  similarly,  the  FE  and  aerodynamic  models 
corresponding  to  the  scaled  structure  will  also  be  different.  The  wind  tunnel  chosen  to  perform  the  tests 
will  constrain  the  dimensions  of  the  model  and  the  speed  it  is  tested  at. 

The  differences  between  the  full  scale  (subscript  s)  and  model  (subscript  m  -  physical  and  computational) 
mean  that  it  is  not  possible  to  simply  scale  the  FE  matrices.  Instead  comparison  must  be  made  of 
quantities  such  as  modal  quantities  (frequencies  and  mode  shapes),  influence  coefficients  and  non- 
dimensionalised  deflections. 

Here  it  is  proposed  to  use  the  following  scaled  quantities. 

1 .  Geometric  Scaling 
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2.  Velocity  Scaling 

y  =n  v 

s  v  m 

(28) 

3.  Air  Density  Scaling 

Ps=npPm 

(29) 

Note  that  the  Equivalent  Air  Speed  could  be  used  in  order  to  eliminate  the  need  for  density  information. 

4.  Natural  Frequencies 

Rewriting  the  aeroelastic  equations  in  non-dimensional  form  involves  the  transformation  into  non- 
dimensional  time.  The  natural  frequencies  of  the  full-scale  and  model  structure  at  the  reference  test 
conditions  must  be  the  same,  thus 
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i.e.  the  reduced  frequencies  must  remain  the  same. 


(30) 


5.  Mode  Shapes 

The  mode  shapes  corresponding  to  the  various  natural  frequencies  of  the  full  scale  and  model  structure 
must  be  the  same.  These  shapes  will  be  complex  to  some  degree  depending  upon  the  characteristics  of 
the  unsteady  aerodynamics  terms. 


6.  Static  Aeroelastic  Deflections 


The  non-dimensionalised  static  aeroelastic  equations  (26)  show  that  equivalent  non-dimensionalised 
deflections  at  the  reference  test  conditions  must  be  the  same.  e.g.  for  the  wing  tip  at  the  reference 
conditions 
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7.  Gust  Response 

Gust  velocities  will  scale  in  the  same  way  as  the  airspeed.  In  order  to  achieve  an  equivalent  scaled 
response,  the  same  maximum  normalised  deflection  at  some  points  of  the  structure  to  the  scaled  gust 
must  be  obtained. 
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8.  Non-Linear  Static  Aeroelastic  Deflections 


Although  the  FE  model  of  a  non-linear  analysis  is  much  more  complicated  than  the  linear  description 
shown  above,  a  similar  scaling  approach  can  be  used.  In  the  same  way  as  for  the  linear  static  aeroelastic 
deflections,  the  equivalent  non-dimensionalised  non-linear  deflections  at  the  reference  test  conditions 
must  be  the  same.  e.g.  for  the  wing  tip  at  the  reference  conditions 
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and  this  will  enable  any  non-linear  deflection  behaviour  to  be  scaled. 


(33) 


9.  Buckling 

The  critical  speed  at  which  buckling  (linear  or  non-linear)  will  occur  simply  depends  upon  the  velocity 
scaling 
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Given  the  above  geometric,  velocity  and  density  scalings,  an  aeroelastically  scaled  FE  model  can  be 
found  by  optimising  the  dimensions  of  the  internal  structure  and  skin  thickness  (and  possible  extra 
distributed  masses)  to  meet  a  series  of  pre-defined  reduced  frequencies,  mode  shapes,  non- 
dimensionalised  static  aeroelastic  deflections  and,  if  required,  gust  and  buckling  behaviour.  It  is  unlikely 
that  an  exact  match  will  be  made  for  all  conditions. 


7.  Example  of  Optimisation  Approach 

The  optimisation  approach  was  developed  using  a  NASTRAN  finite-element  model  of  a  Convair  990  (see 
figure  21).  The  objective  was  to  create  a  2-spar  and  multi-ribbed  wing  that  would  have  similar  (within  5- 
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10%  for  the  first  few  modes)  modal  properties  to  the  original  aircraft  wing,  whilst  not  altering  the  geometric 
plan-form  properties  of  the  aircraft. 


Figure  21 .  Aircraft  FE  Model 


7.1.  Method 

The  original  Convair  990  FE  wing  (see  figure  22)  is  a  representation  of  the  true  wing,  with  the  PBEAM 
elements  (yellow,  running  from  root  to  tip  forming  a  single  line)  having  no  physical,  inertial  or  stiffness 
resemblance  to  the  true  spar.  The  first  step  was  to  reduce  the  FE  model  down  to  the  starboard  wing  so 
that  the  mode  shapes  could  easily  be  identified.  Next,  a  modal  analysis  of  this  wing  was  performed,  the 
resulting  modes  being  those  that  were  to  be  re-created  in  the  modified  model.  The  original  spanwise 
beam  was  then  removed  and  two  new  spanwise  beams  were  added  at  quarter  and  three-quarter  chord 
positions  to  represent  the  two  spars.  Additionally,  6  beams  were  added  to  model  the  ribs  (see  figure  23).  A 
MATLAB  program  was  then  written  to  allow  Genetic  Algorithm  (G.A.)  optimisation  via  a  macro  program 
with  NASTRAN.  Each  spar  was  divided  into  4  property  sets  and  each  rib  formed  a  property  set  resulting  in 
14  property  sets.  The  mass  and  stiffness  properties  of  each  spar  element  and  rib  element  (both 
rectangular  cross-sections)  were  controlled  via  the  variation  of  the  cross-sectional  dimensions  (width  and 
depth),  giving  a  total  of  28  variables  in  the  model.  These  variables  were  restricted  within  a  specified  range 
in  order  to  allow  any  solution  to  be  realistic.  The  G.A.  program  creates  an  initial  set  of  these  variables, 
then  NASTRAN  runs  the  analysis  with  these  values.  The  G.A.  program  then  compares  the  new  modal 
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frequencies  to  the  target  frequencies  before  creating  a  new  set  of  variables  biased  towards  a  better 
solution.  This  process  continues  until  suitable  convergence  is  achieved. 


Figure  22.  Aircraft  Wing  FE  Model 


Figure  23.  Equivalent  Two  Spar  Wing  FE  Model 


7.2.  Results 

Two  cases  were  analysed;  firstly,  using  the  first  three  modes  as  the  optimisation  objective,  and  secondly 
trying  to  re-create  the  first  five  modes  of  the  original  wing.  The  results  are  shown  below  alongside  the 
target  frequencies. 
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Case  1  -  3  modes 


Target  frequencies: 

1 .6829 

3.3127 

5.4791  [Hz] 

Converged  frequencies: 

1 .7893 

3.1507 

5.5651  [Hz] 

Percentage  difference: 

4.9411% 

Case  2  -  5  modes 

Target  frequencies: 

1 .6829 

3.3127 

5.4791  6.0044  7.541  [Hz] 

Converged  frequency:  1 .5705 

3.6555 

4.7721 

6.6425  8.7635[Hz] 

Percentage  difference: 

10.387% 

For  the  first  case  it  can  be  seen  that  the  technique  works  very  well  at  achieving  the  desired  modal 
properties;  the  difference  between  the  results  and  the  original  wing  is  less  than  5%.  As  we  try  to  optimise 
for  a  greater  number  of  modes  (five),  the  analysis  results  in  a  less  accurate  but  still  respectable  similarity 
to  the  original  modes  (approximately  10%  difference). 


8.  Initial  Scaled  Model  Design 

The  above  optimisation  approach  is  to  be  used  for  the  design  of  a  half  sensorcraft  wind  tunnel  model. 
This  will  be  used  by  Lt.Col  Vanessa  Bond  for  her  studies  into  the  non-linear  deflection  /  buckling 
behaviour  and  at  Manchester  as  a  bigger  demonstrator  of  the  gust  load  alleviation  device. 

The  full-scale  sensorcraft  is  shown  in  figure  24  and  the  internal  construction  of  the  proposed  wind  tunnel 
model  in  figure  25.  Note  that  an  outer  metal  skin  will  be  applied  as  well.  The  design  parameters  in  the 
optimisation  process  will  be  the  spar  and  rib  dimensions  and  also  the  skin  thickness.  It  is  worth  noting 
that  the  model  wing  is  attached  in  all  6  DOF  at  the  wall  whereas  the  rear  wall  attachment  allows  some 
vertical  motion  via  a  spring  attachment.  The  two  spar  and  rib  construction  will  be  used  as  it  is  relatively 
simple  to  manufacture  and  through  changing  the  parameters  it  should  be  possible  to  achieve  a  scaled 
model. 
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One  difficulty  with  implementing  the  optimisation  of  the  scaling  process  is  that  the  University  of 
Manchester  is  not  allowed  to  have  access  to  the  full-size  FE  model  and  consequently  have  to  rely  upon 
Vanessa  Bond  at  AFIT  to  perform  the  required  calculations  for  us. 


Figure  24.  Full-size  Sensorcraft  Model 
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Figure  25.  Two  Spar  Wing  Tunnel  Model  Internal  Structure  FE  Mesh 


9.  Conclusions 

A  prototype  wind  tunnel  model  of  a  passive  gust  load  alleviation  device  has  been  designed,  manufactured 
and  tested.  Initial  results  show  that  the  concept  performs  as  expected.  An  approach  for  the  aeroelastic 
scaling  of  a  half  sensorcraft  model  has  been  defined  and  a  methodology  to  achieve  this  demonstrated 
upon  an  example  wing.  This  approach  will  be  applied  for  the  design  of  a  half  sensorcraft  model  to  be  built 
and  tested  at  Manchester. 
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